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Abstract. We briefly discuss models of energetic particle acceleration by supernova 
shock in active starforming regions at different stages of their evolution. Strong shocks 
may strongly amplify magnetic fields due to cosmic ray driven instabilities. We discuss 
the magnetic field amplification emphasizing the role of the long-wavelength instabilities. 
Supernova shock propagating in the vicinity of a powerful stellar wind in a young stellar 
cluster is argued to increase the maximal CR energies at a given evolution stage of super- 
nova remnant (SNR) and can convert a sizeable fraction of the kinetic energy release into 
energetic particles. 



1. Introduction 

The observed spectra of Galactic cosmic rays 
(CRs) are shaped by two basic processes - 
the acceleration in the sources and the subse- 
quent propagation in cosmic magnetic fields 
and radiation fields. A transition from galac- 
tic to extragalactic cosmic rays is expected 
to occur somewhere between 10'^ eV and 
10'^ eV (e.g. Hillas, 2005; Aharonian et al., 
2011). A preponderance of evidence suggests 
that the particle acceleration mechanism most 
likely responsible is diffusive shock accelera- 
tion (DSA) (e.g., Blandford & Eichler, 1987; 
Jones & Ellison, 1991; Malkov & Drury, 
2001). 

An important question for DSA and CR 
origin has always centered around the maxi- 
mum particle energy a given shock can pro- 
duce. For a shock of a given size, age, and 
magnetic field geometry, the maximum CR 
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energy depends mainly on the power in the 
longest wavelength turbulence. The weakly 
anisotropic distribution of accelerated parti- 
cles, i.e., CRs is considered in § 2 as an agent 
producing this turbulence in a symbiotic re- 
lationship where the magnetic turbulence re- 
quired to accelerate the CRs is created by the 
accelerated CRs themselves. 

Other important issue of the high energy 
CR acceleration models is the environment 
where the supernova exploded. Some models 
of CRs acceleration by SNRs in active star 
forming regions are briefly discussed in § 3. 
OB-associations and young globular clusters 
are observed both in the Milky Way and LMC. 
The study of the stellar content of the galac- 
tic object Cygnus OB2 by Knodlseder (2000) 
have revealed that the number of OB mem- 
ber stars can be estimated as large as 2600 + 
400, while the number of O stars amounts to 
120 + 20. The high number of stellar X-ray 
sources detected with Chandra by Wright & 
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Drake (2009) confirmed the status of Cygnus 
OB2 as one of the most massive SFRs in the 
Milky Way. Given the apparently compact size 
of Cygnus OB2 one may expect a number of 
massive stars with strong winds to be in a close 
proximity with less than 10 pc separation in 
addition to the well known coUiding-wind bi- 
naries that are expected to be particle accel- 
erators (e.g. Eichler & Usov, 1993). Another 
very compact young stellar cluster Westerlund 
2 containing more than dozen O stars was 
found in the Carina region with estimated age 
to be younger than 4 Myrs, so the most mas- 
sive stars are expected to explode there within 
the next few Myrs. OB associations with su- 
pernova explosions are creating superbubbles 
(SBs). Bamba et al. (2004) discovered both 
thermal and nonthermal X-rays from the shells 
of the SB 30 Dor C in the LMC. The X-ray 
morphology was reported as a nearly circular 
shell with a radius of 40 pc, which is bright 
on the northern and western sides. Maddox 
et al. (2009) analyzed Suzaku observations of 
the SB around the OB association LH9 in 
the HII complex N 1 1 in the Large Magellanic 
Cloud. Their X-ray spectral analysis revealed 
that the hard X-ray emission (> 2 keV) re- 
quires a hard nonthermal power-law compo- 
nent. The energy budget analysis for Nil us- 
ing the known stellar content of LH9 indicated 
that the observed thermal and kinetic energy 
in the SB is only half of the expected mechani- 
cal energy injected by stars, consistent with the 
expectation of SB models with efficient CR ac- 
celeration (e.g. Bykov, 2001; Butt & Bykov, 

2008) . Diffuse X-ray emission was detected 
from many sites of massive star formation: the 
Carina Nebula, M17, 30 Doradus, NGC 3576, 
NGC 3603, and others (e.g. Townsley et al., 
2011). The H.E.S.S. telescope detected high- 
energy gamma rays from starburst galaxy NGC 
253 supporting the ideas of efficient CR accel- 
eration in active starforming regions (Acero, 

2009) . 

2. Magnetic field amplification in CR 
acceleration sources 

Fast and efficient CR acceleration by Fermi 
mechanism requires ffiat particles are multi- 



ply scattered by magnetic fluctuations in the 
acceleration source (e.g. shock). The magni- 
tude of the required magnetic fluctuations is 
substantially higher than the ambient magnetic 
turbulence forcing a bootstrap scenario where 
the accelerated particles amplify the turbu- 
lence required for their acceleration. The study 
of turbulence generation associated with CRs 
and DSA has a long history. Magnetic field 
amplification due to the resonant cosmic-ray 
streaming instability was studied in the con- 
text of galactic cosmic -ray origin and prop- 
agation since the 1960s (see e.g. Kulsrud & 
Cesarsky, 1971; Wentzel, 1974; Achterberg, 
1981; Berezinskii et al., 1990; Zweibel, 2003). 
It was proposed by Bell (1978) as a source of 
magnetic turbulence in the test particle DSA 
scenario, and nonlinear models of DSA in- 
cluding CR-driven instabilities and magnetic 
field amplification were investigated by Bell 
(2004); Amato & Blasi (2006); ZirakashviH & 
Ptuskin (2008),Vladimirov et al. (2008, 2009) 
and Reville et al. (2009). 

It is instructive to summarize the growth 
rates for magnetic instabilities that the quasi- 
linear theory predicts for weakly anisotropic 
CR distributions of the form 
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where - particle pitch-angle, // - cosO , 
jSj - u jc, with the flow (shock) velocity u,, 
6{p) - is the magnitude of the second har- 
monic anisotropy. We included here the sec- 
ond harmonic anisotropy to illustrate the effect 
of the CR-firehose instability on the magnetic 
field amplification. The substantial anisotropy 
of the second harmonics can be expected in 
more complex plasma flows than just an iso- 
lated shock. The isotropic part of the distribu- 
tion function is a power-law of the form 



Nip) 



Tier {a - 3) 

[1 - A"-3]p' 



(0-3) 



(2) 



where the CR momenta po < p < p„„ A - j^, 
a - is the CR power law index, and «(,,- is the 
CR number density. The CR mean free path 
A = T]rg{p), where the particle gyroradius is 
rgip) - and 77 > 1. Then, following the 
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Fig. 1. The two panels show the growth rates for two opposite circular polarization modes (shown in left 
and right panels) of the parallel propagating modes as a function of the wavenumber. The CR resonant 
streaming instability and the Bell instability The dot-dashed curves effect of short-scale turbulence on the 
long-wavelength instability. The model parameters are A = 100, ySj = 0.01 and a = 4.0. For comparison, the 
dot-dashed curves are calculated for the model of the ponderomotive CR current instability in the presence 
of the short-scale turbulence (Bykov et al., 2011). The dotted curves illustrate the maximal case of the 
CR-firehose instability (i.e., i5 oc /?,). The dashed curve corresponds to (5 = 5/3^. 



standard linear analysis of the kinetic equation 
in the intermediate regime ?; ' < xq < 1 (see 
e.g. Bykov et al., 2011) one may get the fol- 
lowing dispersion relation: 
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where kq = , x - krJp) , xq - 

c Bo 

krgipo) , x„ - krgipin). The signs + corre- 
spond to the two opposite circularly polarized 
modes under investigation. The second term 
in the right hand site in corresponds to the 



fast growing mode discovered by Bell (2004) 
for Xq > 1. The instability results in ampli- 
fication of a mode with wavenumbers ko > 

k > — . For AO < 1 the term in Eq. 3 is re- 

sponsible for the well known resonant stream- 
ing instability (e.g. Achterberg, 1981; Zweibel, 
2003; Pelletier et al., 2006; Marcowith et al., 
2006; Amato & Blasi, 2009, and the references 
therein). The last term in the rh.s. of Eq. 3 de- 
scribed the CR firehose instability. In the long- 
wavelength regime x,„ <k 1 a simplified form 
of Eq. 3 can be derived 



OJ 



1 + 



(JyS^'lnA 
(1-A-') 



(7) 



The growth rates of the unstable modes from 
the dispersion relation Eq. 3 are shown in 
Figure 1. The growth rates of the resonant 
and the Bell instabilities are shown by solid 
lines in Figure 1. Since x„, oc ^ at x„, ~ 1, 
the growth rate Im{io) oc k^l^ - the long- 
wavelength regime pointed out by Schure (& 
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Bell (2011) (for A = 100). At longer wave- 
lengths however the second term in the square 
brackets in Eq. (7) dominate and the CR- 
firehose type instability appear at x,„ << 1 . The 
CR-firehose instability rate is indicated by the 
dashed line at the left panel in Figure 1. The 
CR-firehose instability may even dominate the 
regime Im{<jS) oc k^^'^ in the case of the second 
harmonic comparable with the first harmonic 
(CR current) of the anisotropy, i.e. 6 oc p^. 
The CR-firehose growth rate in that regime is 
shown by the dotted curve in Figure 1. That 
likely requires a more complicated structure of 
the magnetic field in the flow. The flow in the 
vicinity of the SNR colliding with a powerful 
stellar wind (see e.g. Velazquez et al., 2003) 
discussed below as the PeV CR acceleration 
site may be relevant. To derive the growth rates 
of the modes in the long-wavelength regime 
^A < 1 (i.e. .^o < ?7 ') the collisionless kinetic 
equation approach discussed above is not ap- 
propriate. One should average the kinetic equa- 
tion for the relativistic particles, the equations 
of the bulk plasma motions and the induction 
equation over the ensemble of the short scale 
fluctuations produced by CR instabilities in the 
collisionless regime e.g. by the fast Bell insta- 
bility. In the presence of the short-scale fluc- 
tuations, the momentum exchange between the 
CRs and the flow in the hydrodynamic regime, 
results in a ponderomotive force that depends 
on the CR current in the mean-field momen- 
tum equation of bulk plasma (Bykov et al., 
201 1). As a result, there exist transverse grow- 
ing modes with wavevectors along the initial 
magnetic field with growth rates that are pro- 
portional to the turbulent coefficients deter- 
mined by the short scale fluctuations. The mag- 
netic field amplification in that regime only 
weakly depends on the shock velocity [as it fol- 
lows from Eq. (48) in Bykov et al. (201 1)], that 
is important for the evolution of the maximal 
energy of CRs accelerated by DSA. 

3. Particle acceleration by a 
supernova remnant in OB 
association 

The maximal energies of CRs accelerated by 
supernova shocks with magnetic field amplifi- 
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Fig. 2. Simplified view of the expanding supernova 
remnant approaching a powerful stellar wind. The 
configuration is shown to be favorable for 10'^ eV 
regime CR proton acceleration at a few thousands 
years old SNR in SFRs. i?„. and R^,, are the radii of 
the stellar wind termination shock and SNR forward 
shock, respectively. k„, is velocity of the stellar wind, 
Usn is the SNR forward shock velocity, L is the dis- 
tance between the two flows. 
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Fig. 3. CR particle spectrum (solid line) from the 
time-dependent DSA simulation of the SNR for- 
ward shock approaching the O-star wind termina- 
tion shock [see Figure 2]. For the comparison the 
CR spectrum at solitary shock is shown as the dot- 
ted line curve. 



cation can reach the knee region of the GCR 
spectrum (e.g. Blandford & Eichler, 1987; 
Bell, 2004; Berezhko & Volk, 2007; Caprioli 
et al., 2011). Moreover, it was recently argued 
by Ptuskin et al. (2010) that the maximum en- 
ergy of accelerated particles differs strongly for 
different supernova types. The energy may ex- 
ceed 10'^ eV for protons accelerated in an iso- 
lated Type lib SNR. The highest energies CRs 
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are thought to be accelerated before the Sedov 
stage in isolated SNRs. 

For supernova exploded in a young stellar 
cluster the presence of a nearby strong stellar 
wind allow to increase the maximum energy 
of accelerated particles comparing to that in an 
isolated SNR. The SNRs can produce the high 
energy CRs at the Sedov evolution stage. 

The most massive O stars begin to explode 
as core collapsed supernovae about a few mil- 
lion years after the formation of a young stellar 
cluster. The young stars with masses above ~ 
16 Mq (of BO V type and earlier) are thought 
to create hot low-density bubbles and HII re- 
gions with radii ~ lOpc surrounded by a mas- 
sive shell of matter swept up from the parent 
cloud by the stellar wind. In this case, a strong 
supernova shock propagates for a few thousand 
years in tenuous circumstellar matter with a ve- 
locity well above 10^ km/s before reaching the 
dense massive shell. The blast wave of the SNR 
is expected to accelerate the wind material pro- 
ducing ultra-relativistic ions and electrons. The 
accelerated particles escape from the SNR for- 
ward shock presumably at the highest energy 
regime (e.g. Ptuskin et al., 2010; Gabici, 201 1; 
Ellison & Bykov, 201 1) and can reach the ter- 
mination shock of the stellar wind of a nearby 
massive star. We modeled the energetic particle 
acceleration in the region where the expanding 
supernova shell is approaching a powerful stel- 
lar wind of a young massive star as it is illus- 
trated in Figure 2. At the evolution stage where 
the mean free path of the highest energy CR 
is comparable to the distance between the two 
shocks L the system is characterized by an un- 
usually hard spectral energy distribution illus- 
trated by the solid line in Figure 3. The stage 
can last more than 1,000 yrs. For comparison 
we showed with the dotted line the CR spec- 
trum accelerated by an isolated SNR shock of 
the same age. The case of Bohm diffusion was 
simulated with D{p) oc p. 

A simple analytical kinetic model of par- 
ticle acceleration by the Fermi mechanism in 
the converging flows carrying fluctuating mag- 
netic fields amplified by anisotropic CR distri- 
butions can be considered. It is a model of high 
energy CR acceleration by a SNR expanding 
in a compact OB-association. At some expan- 



sion phase the distance L between the SNR ap- 
proaching the stellar wind shock is less than the 
mean free path of the highest energy CR parti- 
cle in the SNR shock precursor Then CR dis- 
tribution function around the shocks (/ = 1 , 2) 
can be approximated as 

fi {x, p, t) = Ap-^ exp {-%■ |x|) X 

XH(p-pQ)H(t-Ta) 

where the highest energy CR acceleration time 




Po 



The acceleration mechanism can provide 
efficient creation of a nonthermal particle pop- 
ulation with a very hard energy spectrum, con- 
taining a substantial part of the kinetic energy 
released by the supernova. The high energy 
nonthermal emission of the sources is charac- 
terized by a very hard spectral energy distri- 
bution peaked at the maximal photon energies 
and have the apparent properties of the "dark 
accelerators". The TeV source in the vicinity 
of the y-Cygni SNR can be belong to the class. 

At the later stage of the young stellar clus- 
ter evolution multiple supernova explosions 
with great energy release in the form of shock 
waves inside the superbubbles are argued as a 
favorable site of nonthermal particle accelera- 
tion. The collective effect of both stellar winds 
of massive stars and core collapsed supernovae 
as particle acceleration agents were discussed 
by Bykov & Toptygin (2001); Parizot et al. 
(2004); Torres & Domingo-Santamaria (2007); 
Ferrand & Marcowith (2010). 
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